function ZCRanalysis(noWords,nRep,frameSize,overlap);

%Autor: Paula Sanz Leon
%Fecha: Diciembre 2008

%Esta rutina calcula la media de ZCR para poder estimar el umbral para el
%algoritmo de VAD.

for	i = 12:noWords
    for j = 1:nRep
    
        q           = ['/home/paula/espich/TrainingData/Speaker/spkr_' num2str(i) '_' num2str(j) '.wav'];   %DIRECTORIO DONDE SE ENCUENTRA LA SECUENCIA DE ENTRENAMIENTO
        [signal,fs] = wavread(q);
        signal      = signal / max(abs(signal));     %normalizacion
        L           = length(signal);
        frameMat    =buffer(signal, frameSize, overlap);
        zcr1        =sum(frameMat(1:end-1, :).*frameMat(2:end, :)<0); %Método 1
        zcr2        = sum(frameMat(1:end-1, :).*frameMat(2:end, :)<=0);%Método 2, incluye los ceros
        time        = (1:L)/fs;
        frameNum    = size(frameMat, 2);
        frameTime   = ((0:frameNum-1)*(frameSize-overlap)+0.5*frameSize)/fs;
        figure;
        subplot(3,1,1); 
                        plot(time, signal); 
                        ylabel('Amplitud Normalizada', 'FontSize', 12);
                        xlabel('Tiempo(seg)', 'FontSize', 12);
                        title(['ZCR', sprintf('\n'),q], 'FontSize', 12);
        subplot(3,1,2); 
                        plot(frameTime, zcr1, '.-', 'Color', [0, 0.75, 0.75]);
                        ylabel('Cruces por Cero', 'FontSize', 12)
                        title('Metodo 1', 'FontSize', 12);
        subplot(3,1,3);
                        plot(frameTime,zcr2,'.-','Color',[1, 0.2, 0.4])
                        ylabel('Cruces por Cero', 'FontSize',12);
                        xlabel('Tiempo(seg)', 'FontSize',12);
                        title('Metodo 2', 'FontSize', 12);
        
        ZCR1(j,i)     = mean(zcr1);
        ZCR2(j,i)     = mean(zcr2);
    end
     
    AVzcr2(i)   = mean(ZCR2(:,i))
    AVzcr1(i)   = mean(ZCR1(:,i))
end

